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We propose a one-dimensional Hamiltonian H\d which supports Majorana fermions when d^2_y2- 
wave superfluid appears in the ultracold atomic system and obtain the phase-separation diagrams 
both for the time-reversal-invariant case and time-reversal-symmetry-breaking case. From the phase- 
separation diagrams, we find that the single Majorana fermions exist in the topological superfiuid 
region, and we can reach this region by tuning the chemical potential jj, and spin-orbit coupling 
a_R. Importantly, the spin-orbit coupling has realized in ultracold atoms by the recent experimen- 
tal achievement of synthetic gauge field, therefore, our one-dimensional ultra-cold atomic system 
described by Hm is a promising platform to find the mysterious Majorana fermions. 
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1. INTRODUCTION 

Since the discovery of the fractional quantum hall state 
[1], the concept of the topological order, which was first 
proposed explicitly by Wen [2], has been developed very 
fast and used in many condensed matter systems. A 
gapped system, such as the Pfaffian state proposed by 
Moore and Read [.3], which possesses topological order, 
may have practical use in the topological quantum com- 
putation (TQC) on account of its quasi-particles' non- 
trivial properties, such as non-Abelian statistics, and the 
tolerant ability to the decoherence from the environment 

Read and Green [5] pointed out that the zero energy 
Majorana fermion modes exsiting at the cores of a 2D 
spinless p-wave superconductor in the weak-phase [6, 7] 
are the same as the nonabelions in the Pfaffian state 
[3], and they are non-Abelian quasi-particles [8]. Ki- 
taev constructed a toy model and showed that Majorana 
fermions exist as end states of a spin-polarized 1-D super- 
conductor [9]. This model supplies an insightful way to 
find the interesting single Majorana fermion. Recently, 
many groups have proposed different systems to engi- 
neer topological superconductivity (TSC) with Majorana 
fermions as a bonus [10-17]. Among them, papers [14, 15] 
recognize that the topological superconductivity can be 
perhaps most easily engineered in one-dimensional semi- 
conducting wires deposited on an s-wave superconduc- 
tor, and provide the first realistic experimental setting 
for Kitacv's model and a platform to find and manipu- 
late single Majorana fermion by braiding [18]. The au- 
thors of [in] propose that Au wires in proximity to doped 
ljSCO{L2-xSrxCu04), a ((i^2_y2) wave superconductor 
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can be a more promising candidates for realizing single 
Majorana fermion. 

In addition to fractional quantum hall systems and 
topological superconductors, topological non-trivial su- 
perfluid, which can be deduced from an underlying nor- 
mal superfluid, e.g., s-wave superfluid, also supports non- 
Abelian Majorana fermions. With the rapidly developing 
technology available for the quantum control of ultra-cold 
atomic systems, clean environment and highly tunable 
parameters, ultracold atomic systems may serve as an 
idea platform for the observation of topological superflu- 
idity and topological phase transition. Importantly, the 
realization of spin-orbit coupling in ultracold atoms by 
the recent experimental achievement of synthetic gauge 
fleld has made a flrm step to engineer topological super- 
fluidity and non-Abelian quasi-particles therein [20, 21]. 

With the introduction of the spin-orbit coupling, the 
energy gap of superconductors or superfiuids with asym- 
metry pairing wave functions will close at some points in 
the first Brillouin zone. If the system is under certain 
symmetries, e.g. time reversal symmetry, particle-hole 
symmetry etc., and the manifold is not closed, the system 
will possess robust gapless edge excitations protected by 
these symmetries. However, once the symmetry protect- 
ing the gapless excitations is broken, e.g. Dili class to D 
in 3D case [22], the gapless excitations will no longer be 
protected and will be absorbed by some random impuri- 
ties. Luckily, there are some cases [22] even the symmetry 
is broken, there will be gapless excitations that can still 
exist robustly. 

In this article, we study the one-dimensional ultra- 
cold atomic system with spin-orbit coupling both for the 
time reversal symmetry and the time-reversal-symmetry- 
breaking case. We find that there is only one topologi- 
cally nontrivial superfluid phase, which is always pro- 
tected by an energy gap away from the normal superfluid, 
which agrees with the result in Ref.[19]. What most in- 
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tcrests us is that we can reach the topological superfluid 
(TSF) region explicitly by tuning the parameters accord- 
ing to the phase-separation diagrams and directly probe 
the single Majorana fermions. 



2. MODEL STUDY 

We consider a one-dimensional ultracold atomic sys- 
tem with spin-orbit coupling and the Hamiltonian is 



Hid — Ht + Hso + Hj + Hz 



(1) 



where 



j.a j,a 

Hso ^-\Y1 (^"fl^]+i,a(c^y)a/3V'j,/J + h.c. 

Q/3 ij 

Hz^J2^z{4'U^,,t-^li^j,l 



(2) 



where ipj is a fermion operator at site j, a and /3 are the 
spin indices, t is the hopping amplitude, fi is the chemical 
potential, an is the spin-orbit coupling strength, Hj is 
the interaction, gij is the interaction strength between 
site i and site j, and is negative in this model, ay is a 
pauli matrix, Hz is the Zeeman term which breaks time 
reversal symmetry, and Vz denotes the strength. 

By Fourier transformation and mean field approach, 
the above Hamiltonian will take the form 



Hu 



Ht + H. 



so 



H. 



SF 



H, 



(3) 



where 



Ht = 
Hso 



^[tC0s(fc)-t-^] -0j „V'fc,a, 

- iaii sin(/s) (V'l.tV'fe,; - ^l,i^kA^ : 

k 

HsF = J2 [Mk)^U^U^ 

k 

-HA*(fc)V'-fc,i^fe,t]-|AolV>/, 

Hz = Y.^z(i^Ui^k.^-^i^^k.i), (4) 



where A(fc) = j^Y.k' dik ~ k' ){%!:_ ^> = Aq cos(fc) 

is the pairing amplitude, g{k — k ) = gcos{k — fc ) is the 
Fourier form of gij, g is a real constant. Aq is a complex 
constant, here we make it real for convenience. Here the 
pairing we are interested in is the d^^-y^ type. N is the 
number of sites and the lattice constant a is set as unit, 
J = . Strictly speaking, making mean field approx- 
imation is not proper for a one-dimensional system, as 
fiuctuations are strong and there is no true long range 
order in a homogeneous system with non-zero tempera- 
ture in the thermodynamic limit, according to the well- 
known Hohenberg-Mermin- Wagner theorem. However, 
as shown in Ref. [2.'!], by confining the system in a box 
with finite length L or in a harmonic trap to avoid this 
technology difficulty, the authors found the mean-field 
methods provide a useful description in the weakly or 
moderately interaction regimes by comparing the mean- 
filed result with the exact of asymptotically exact Bethe 
Ansatz solutions. In the following, we will set the length 
of the system to be L = Na = N and T = 

In the momentum space, under the Nambu spinor 
representation ^'(fc)t = {V'fc .|., -^l^.j., V'-fc.j., V'fc.tli the 
Hamiltonian can be rewritten as 



H,D = \Y.'^{k)^H{k)-^{k)~^l/J, 



(5) 



where 



H{k) 
h{k) 

m 



h{k) A(fc) 
A(fc)t -h^{-k) 

Ekcro + ctR sin(fc)(Ty + Vz(Jz 

iA{k)<7y. 



(6) 



where Sk — —tcos{k) — fi. After diagonalizing, the Hamil- 
tonian is of the form 



= I E {[(^i(^) - ^3(fc)] aI,t«fe,T + [E2ik) - E^{k)] + E^[k) + i?4(fc)} 



A^ 
J 



(7) 
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where Ei{k), E2{k), E3{k), E4{k) are in the order {++,H — , — h, } of 



E{k) ^ ±Jel + al^ik) + A2 cos2(fc) + Vi ± 2^ elals\n\k) + e^yj + y|A2 cos2(fc), (8) 



and elhpsis stands for the terms which are constant and 
independent of Aq, a\ -^{Pl ^) is the creation operators of 
the excitation. The ground state is | 0), which satisfies 
ckfc,t(/5fe4) I 0) = Oj the energy of the ground state is 



E, = \Y.^E:,{k)+Em-^, 



(9) 



which needs to satisfy the condition of mcan-field- 
approximation assumption 



2 _ 1 ^ 



A- 



1^2 



E-i(k)A Ei{k)A 



cos^{k), 



(10) 



where A = J sla% sin'' {k) + ejVj + VjAg cos2(fc). 



3. RESULTS AND DISCUSSIONS 

First, for the time- reversal invariant case (Vz = 0), we 
have 



E3{k),Ei{k) - -y^(£fc ± an sin(fc))2 + Ag cos2(fc). (11) 

When Vz = 0, the Hamihonian H{k) possesses a chiral 
symmetry. Therefore, the Hamihonian can be brought 
into an off-diagonal form by a unitary transformation 

[24] 



Ql 



Qk 



(12) 



with 



1 



1 



^/2 W 1 M -^T 



(13) 



and T = lay, Qk = i[e^«-('=H'^o - ^y) + e''+(^)(ao + 

a )\ where e'®*''"") = -tcos{k)-^l±a^tsin{k)+i^locos{k) 

' y/[-t cos(fe)-M±aH sin(fe)]2 + [Ao cos(fe)]2 ' 

In Eq.(12), the meaning of ~ is the magnitude of the 
eigenvalues of H has normalized. 

In one dimension, the Fermi-surface topological invari- 
ant (FSTI) of a time-reversal-invariant (TRI) supercon- 
ductor or supcrfluid is [24] 



NiD = W[sgn{5s 



(14) 



where s is summed over all the Fermi points between 
and TT. If a system has odd number of Fermi points 
between and tt with negative pairing, in other words, 
NiD = — 1, the system is non-trivial, otherwise A^id ~ 1, 
and the system is trivial. 

In our system, there is only one Fermi point {k ~ ^ 
with = an) between and tt, so the sign of 6s directly 
determines the Hamiltonian is trivial or non-trivial. The 
sigh of Ss is positive or negative corresponding to the 
change in 9±{k) across A: ^ is — tt or tt [24]. Here, 
in the weak pairing limit, by increasing k from ^ — e to 
^ -I- e with Qffl g (ImI ~ ImI + (e; ^ sire small positive 
constants), we find, for an > ^ (we only consider ^ > 
on account of symmetry) , the real and imaginary part of 
—t cos(fc) — n + aR sin k + zAq cos(fc) change as 



with 



-(fc) — > 0+{k)-7r 
Ae+{k) = +7r 



-(/c) + 27r, 



for aR < fi, 



6 + ie — > —S — ie — > 6 — ie 



with 



9+{k)^ 0+{k)-2n 
Ae+{k) = 



9+(fc)+7r, 



(15) 

(16) 
(17) 

(18) 



((5, e are small positive real numbers) therefore, when fi < 
c^Ri Nijj = sgn{Ss) = —1, the system is non-trivial; when 
/i > afi, NiD = sgn{5s) = 1, the system is trivial. 

Base on Eq.(lO), Eq.(ll) and Eq.(14) and the above 
analysis, if we fix the strength of spin-orbit coupling, 
we find that there exists a critical chemical potential /ic 
where a superfluid-normal state transition takes place, as 
shown in Fig. 1(a). Above the critical value, it's the nor- 
mal state where the mean-field-approximation assump- 
tion is not valid. Below it, the superfiuid phase appears 
and exists only when the chemical potential is not too 
large. In this region, when > ur^ the Z2 invariant 
NiD ~ 1 where is the normal superfiuid phase region, 
and when /i < an, Nid = —1 where is the topologi- 
cal superfiuid phase region, as shown in Fig. 1(a). The 
line fi = UR separating NSF from TSF is a critical line, 
crossing this line, the topological number changes and the 
topological phase transition takes place. Because there 
is no symmetry breaking while crossing the critical line. 
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the topological order Nuj is the only parameter to dis- 
tinguish these two phases. 

Second, for the time-reversal-symmetry-breaking case 
{Vz ^ 0), the Z2 Majorana number A4 oiHik) is the new 
topological invariant to determine the system is topolog- 
ically trivial or topologically non-trivial. Following the 
Refs.[9, 19, 25] the Majorana number can be obtained by 
the formula 

M = sgn[PfB{Q)]sgn[PfB{'K)] = ±1, (19) 

where ±1 corresponds to topologically trivial and non- 
trivial states and the antisymmetric matrix B{k) is de- 
fined as B{k) = HiD{cx®cro)- In terms of the parameters 
of the Hamiltonian, the Majorana number can be written 
as 

sgn[{t + ^,Y - (y| - Alm-t + ^if (y| - a^)], 

(20) 

so the Majorana number is = —1 when \ \JV^ — A§ — 
t\ < < \ yJVl - + 1| and 7W = 1 otherwise. 

In this case, we find that the phase diagrams will have 
much more structures than the case discussed above. We 
find that the pairing amplitude Ag monotonically de- 
creases with the increasing Vz^ and when the Zeeman 
field is strong enough, Ag will be dramatically suppressed 
and no longer monotonically increase with | J|, the attrac- 
tive interaction strength, as shown in Fig. 1(b). This can 
be explained by the fact that with the Zeeman field in- 
creasing, the polarization becomes stronger and stronger, 
and finally, this will destroy the singlet dj.2_y2 pairing. 
We also find that n and an have the same effects on 
the pairing amplitude as the Zeeman field, as shown in 
Fig.l(c)(d). 

In the following, we will discuss the most interest- 
ing part of our work. According to the phase dia- 
gram Fig. 2 (a) and Eq.(20), there is a critical line (red 
sohd line), corresponding to ^ = \/V^ — Aq -|~ t and 
^i = t- yJVl - Ag, that separates TSF from NSF. Here 
TSF means that single Majorana fermion exists in this 
region, while NSF means superfluid with Majorana dou- 
blets (as Vz < \/t'^ + A§), which are no longer topologi- 
cally protected due to time-reversal symmetry breaking. 
There is also a critical line (black solid line) , correspond- 
ing to Aq = 0, that separates NSF from N. In Fig. 2(a), 
there is a dashed line. This dashed line corresponds to 
II = \/Vz + o?R- On this line, the energy gap closes, 
however, from Fig. 2(a), we see the Majorana numbers 
on both sides of this line are equivalent, either both are 
trivial or both are non-trivial. Therefore, there is no 
topological phase transition across this line. According 
to Fig. 2(a), we can tune iiji to (0.1, 1.45), aB.lt to (0, 1) 
and reach the TSF region, then we can apply the spatially 
resolved radio-frequency spectroscopy, which would show 
a well isolated signal at zero energy, to detect the associ- 
ated single Majorana fermions. 




12 3 12 1.4 1.6 1.8 

n/t |J|/t 




|J|/t |J|/t 

FIG. 1: (Color online) (a) The parameters of H\u are: J = 
— 1.2, Vz ~ 0,t = 1, there are two critical lines, which separate 
the normal superfluid (NSF) from the topological superfluid 
(TSF) and normal state (N) respectively, (b) t — 1, jj. = 
0.45, Qfl = 0.3, changing the Zeeman field Vz from 0.9 to 
1.05, we see Aq monotonically decreases with the increasing 
Vz, and when the Zeeman field is strong enough, Aq will be 
dramatically suppressed and no longer monotonically increase 
with J, the attractive interaction strength, (c) t — l,Vz = 
0.9, OR = 0.3; (d) t = l,Vz = 0.9, fi = 0.45, (c), (d) are quite 
like (b), Ao also monotonically decreases with the increasing 
/i and a_R, and there are also critical value Hc and Qij^, which 
dramatically suppress Aq. 

In order to obtain the information that how the TSF 
region changes with other tunable parameters and make 
single Majorana fermion measurements more achievable, 
Figs.2(b)-(d) are given. Among them. Fig. 2(b) and 
Fig. 2(c) show that by increasing the attractive interac- 
tion or reducing the Zeeman field, although the NSF re- 
gion are broadened, the TSF region are greatly reduced. 
However, by increasing both the attractive interaction 
and the Zeeman field, as shown in Fig. 2(d), the TSF 
region are greatly broadened. Such effect will be quite 
useful in experiments, as the larger the TSF region is, 
the more detectable the single Majorana fermions are. 
From Figs.2(a)(b)(d), we can find TSF tends to form in 
the high polarization area where the pairing amplitude 
is small and the imbalance of the chemical potential is 
large. 

4. CONCLUSIONS 

In this paper, we propose the one-dimensional Hamil- 
tonian Hid and discuss it both for the time-reversal- 
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our one-dimensional ultra-cold atomic system described 
by Hid is a promising platform to find the mysterious 
Majorana fermions. 

This work is supported by NSFC Grant No.l0675108. 
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FIG. 2: (Color online) The parameters of the Hamiltonian are: 
(a) t=l,Vz ^ 0.9, J = -1.2; (b) t = l,Vz = 0.9, J = -1.5; 
(c) t = 1, Vz = 0.6, J = -1.2; (d) t = l,Vz = 1.2, J = -1.6. 
the red solid line is a critical line, separating TSF from NSF, 
the dotted line is a natural extension of the red solid line, it is 
also a phase boundary, separating TSF from NSF; the black 
solid line is also a critical line, separating NSF from N. The 
dashed line is not a critical line, across it, no phase transition 
takes place. Comparing (b)(c) with (a), we find that the 
TSF regions shrink either with Vz decreasing or with | J 
increasing. In (d), we fix the ration Vz/J to be the same as 
(a) and keep t = 1, then we find the TSF region is broadened. 



invariant case and the time-reversal-symmetry-breaking 
case. By numerical solving the self-consistent equation 
Eq.(lO), we obtain different phase-separation diagrams 
under different conditions. From the diagrams, we find, 
with the spin-orbit coupling and the Zeeman field, TSF 
exists. By tuning the parameters, such as /i and an, 
we can reach the TSF region, where single Majorana 
fermions exist. To detect single Majorana fermions, we 
can apply the spatially resolved radio-frequency spec- 
troscopy, which would show a well isolated signal at zero 
energy. 

With the rapidly developing technology available for 
the quantum control and the introduction of spin-orbit 
coupling to ultra-cold atomic systems, we believe that 
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